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Abstract 



We give a general treatment of twist-three effects in two-photon reactions. We address the 
issue of the gauge invariance of the Compton amplitude in generalized Bjorken kinematics and 
relations of twist-three 'transverse' skewed parton distributions to twist-two ones and interaction 
dependent three-particle correlation functions. Finally, we discuss leading order evolution of twist- 
three functions and their impact on the deeply virtual Compton scattering. 
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1 Introduction. 



Exclusive two-photon reactions involving two hadrons (scattering or production), like deeply vir- 
tual Compton scattering (DVCS) [1], |2|, |3|], 7*7 — > hh with small invariant mass of the hadron 
system |TJ, etc., are of special interest since they involve new nonperturbative characteristics 
which generalize conventional parton densities and/or distribution amplitudes. The leading twist 
factorization gives the amplitude of these processes in terms of a perturbatively calculable coeffi- 
cient function and a skewed parton distribution (for DVCS) ||, |J, |J or a generalized distribution 
amplitude (for 7*7 — > hh) J7| which are responsible for soft physics. However, the twist-two anal- 
ysis of the Compton scattering at non-zero t-channel momentum transfer, A 7^ 0, is inadequate 
since it violates explicitly the gauge invariance of the amplitude due to approximation involved, 
Aj_/ v ' —q 2 <C 1. This calls for a consistent treatment of the effects suppressed in Aj_, i.e. twist 
expansion^. As a first step one takes the contributions linear in A^ which are of twist-three. 
Obviously, then the amplitude will be gauge invariant to the twist-four accuracy. Repeating the 
steps one improves the amplitude accordingly. 

Let us demonstrate the peculiarities of the off-forward kinematics on a simple example of a 
free Dirac fermion theory. As we will see, even and odd parity structures 'talk' to each other in 
the case at hand, when the operators with total derivatives within the context of the operators 
product expansion are relevant. This will be the source of restoration of the electromagnetic 
gauge invariance of the two-photon amplitude defined by a chronological product T {j fl (x)j v (y)} 
of currents j^x) = ip^^ip- The leading light-cone singularity (x—y) 2 — > arises from the hand-bag 
diagram and reads 

T {j^(x)j u (y)} = #(x)7 M) £(x - yh^(y) + iip(yhvfi(y - xfr^x), (1) 

where fi(x) = 7^2 is the free quark propagator. Taking into account the equation of motion 
fttp = and (f) fi(x) = —iS(x) it is a simple task to show that the hand-bag diagram respect current 
conservation. After performing the decomposition of the Dirac structure in Eq. ([I]) it reduces to0 

T{3n(x)jv{y)} = S^pJSpix - y) {ip(x)^{y) - (x <-> y)} (2) 
- ^^u P JS p (x - y) {^(x)7 CT 7 5 ^(l/) + (x *-> y)} , 

where Sp U . pa = g w g va + g^Qup — g^uQpa- One finds that current conservation does not separately 
occur in both terms. Employing the density matrix \Pi)(P2\ for free fields, we find, of course, 
that the amplitude (P 2 \T {j^(x)j u (y)} \P%) respects current conservation, however, there appears 

^^Here we imply kinematical definition of twist. 

2 We use the conventions for Dirac and Lorentz tensors from Itzykson and Zuber M . 
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a cancellation between both matrix elements of different parity on the r.h.s. in Eq. (0). So we 
realize that the current conservation is not manifest in the decomposition (fj). This arises due to 
operators with total derivatives in twist decomposition of ^(x)7 M (l, 75)^(2/), see later Eq. (|28|) . 
However, once this decomposition is performed the gauge invariance is restored automatically. As 
compared to this simple example the only complication in QCD to leading order emerges due to 
the presence of the interaction dependent three-particle contributions. This does not present any 
difficulty and will be solved in the next section. 

Our consequent presentation is organized as follows. The next section is devoted to a detailed 
description of the twist-three formalism for the generalized Compton amplitude and as a result 
the restoration of the gauge invariance sketched above. In section |3| we address the issue of a twist 
decomposition of light-ray (and local) operators in a situation when total derivatives are relevant. 
Then we turn in section 4 to the two cases of matrix elements of operators sandwiched between 
spin-0 and | hadrons and derive relations for twist-three two-particle 'transverse' functions in 
terms of twist-two (Wandzura-Wilczek contribution) and three-particle correlation functions. In 
section 5 we comment on the phenomenological consequences for different asymmetries measurable 
in the DVCS process, and then in section || point out that the evolution of twist-three skewed 
distributions is known in leading logarithmic approximation from analogous studies for forward 
kinematics. Finally, we give our conclusions. 



2 Generalized two-photon amplitude. 

The amplitude of scattering of two virtual photons with momenta q\ and q 2 on a hadron target, 
with incoming (outgoing) momentum Pi (P2), is given by a Fourier transform of a correlator of 
electromagnetic currents, 

T^u = % J d 4 xe^ x (P 2 \T {j tl (x/2)U-x/2)} \P 1 ). (3) 

One introduces the vectors q = \{qi + q 2 ) and P = Pi + P 2 , A = P 2 — Pi = qi — q 2 to describe the 
amplitude. They can be used to construct a pair of the light-cone vectors and n*, such that 
n 2 = n* 2 = and n • n* — 1, as follows 

I 97j / Ml! / Qu T / 111 

^1-4(^)2 2<r»ay 1 - 4(£5) 2 v/l-4(^) 2 4 v /l-4(£<5) 2 

(4) 

where 5 2 = (M 2 — A 2 /4)/q 2 . The transverse metric and antisymmetric tensor are defined as, 
g^ v = g^ u — n^nl — n u n*^, e^ u = e fJiU ^ + . Here and in the following we use the conventions for 
the generalized Bjorken variable £ = —q 2 /P ■ q and skewedness 77 = A ■ q/P ■ q which are scaling 
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Figure 1: Two- and three- leg coefficient functions (plus crossed diagrams) for DVCS process 
which form a gauge invariant amplitude to twist-three accuracy. 



variables of (|3]). Neglecting the corrections 0(5 2 ) we have then to twist-four accuracy 



n, 



-1(2^ + ^), 



n* — -P 

'V 2 



(5) 



This approximation will be used throughout in our analysis. 

In the generalized Bjorken kinematics — q 2 — > oo, P ■ q — > oo, A ■ q — > oo, A 2 = finite, £ and 
being fixed the contribution to (|3|) comes from the diagrams in Fig. |] 



(6) 



where T 2 V comes from hand-bag diagram (left) while T^ v from antiquark-gluon-quark one (right). 
Explicit calculation, done in the light cone gauge B + = 0, gives 



T lv = / dx \ & + \ x >€) \ S wp- Vk p( x >v) + 2ie flupa q a A O p (x,r]) 



2q 2 



for the first one, and 



-Sp V - p+ v O p (x, if - ie pup - [2x A O p (x, if - A K P 



X,T] 



T, 



+ C ( -\x,£f)Spp. v - v S~(x,x -t,t}) - v S+(x + T,x,ri) 



(7) 



for the second. Here S pU]p(7 tensor defined in the introduction and tree coefficient functions are 



C (± \x, = (1 - */t ~ My 1 ± (1 + x/i - iO)- 1 . 



(9) 



We use here the following conventions for generalized functions 



r O p (x, V ) = J ge*" (P 2 |^ p (f,-f) \P 1 ), T K p (x, V ) = J ^e*" (P a |%„(f,-f) [P,), 



I Sf(Xi,X 2 ,T]) 



2tt 2tt 



e2' 



(10) 
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where the two-quark operators are defined by 

t O p (k, —k) = , ip(—Kn) I r p ; ip(Kn), T JC p (k, —k) = ^{—nn) J r + i d p ip(nn), 

with T r = {7 P , 7 P 75} for I = V,A and antiquark-gluon-quark ones read 

v Sp(Ki,K 2 ,K 3 ) = ig4>(n 3 n) j + G +p (K 2 n) ±ij+j 5 G +p (K 2 n) 4>(nin), 
a 5^(«i, «2, k 3 ) = ig${K 3 n) [7 + 7 5 G +p (K 2 n) ± i7 + G +p (/« 2 n)] 

The latter two are related to each other by the 'duality' equation 

ieL*& = ± v S£. 



11) 



per a 



(12) 



(13) 



By means of the quark equation of motion, T^ip = 0, we can obtain the following relation 
between the correlation functions introduced so far 

^ v Oj (k, -k) - i(+d + A 0$ («, -«) + z% («, -«) + ze£ r # A 0+ («, -«) 



+ / dA {w(A - k) v S~ (k, A, -k) + tu(A + («, A, -k)} = 0, 



5+^ («, -k) - («, -«) + iX («, -*) - (*, -*) 



+ / d\ [w(\ - k) 1 ^" (k, A, -k) - w(A + (k, A, -k)} = 0, (14) 

where iw(ac) = —0(k) for the ML prescription on the infrared pole in the gluon propagator. Here 
d stands for the total derivative. In terms of the momentum fraction space functions introduced 
in Eqs. ( |10D these relations read 

2x v O p L (x, V ) + 2r ] ief (T A 0^x,r ] ) - ie^A^O+ix^) - 

dr 



J T -iQ { VS p ( x > x ~ T ,v) + V Sp(x+r,x,r])} = 0, 



r-iO 

2xiej a A Oj(x, ry) + 277 v Oj(x, 77) - v O + (x, 77) - ie^K^x, 77) 

~~Q { VS p( X > X ~ T '^) ~~ ^pfa + r ' X '^)} = °- ( 15 ) 



+ 



So that expressing If in terms of other correlation functions in Eq. fl7|) we get the result 
T, v = -- 2 J dxlrff C^(x^) V + (x,v) + T^ cM(x,t) A + {x,rj) 



+?$ P C { ~\x } V Of(x, 77) + 7^ C«(s, 77) 



(16) 



where the support properties of distributions restrict the integration range of the variable x within 
the interval [—1,1]. Here the Lorentz tensors are 



T$ P = + (2^ - i)K ) g w + (q, + (2f + 77)71*) g up , 

= *W<^ _ 0™ + < ^m) • ( 17 ) 

These expressions are obviously target independent. In the case of scalar target our result reduces 
to the one obtained in Ref. ||. 

It is an easy exercise to check that, after accounting for twist-three corrections, gauge invariance 
is fulfilled up to C(Aj_), i.e. twist-four effects: q\ v T$ = q2pTjg) = 0(A\), where we have used 
Sudakov decomposition for the photon momenta = — ^n p — (£ — r/)n* + ~A^ and q 2p = 
~2ti n P ~ (£ J ^ 1 l) n ]i ~ Interesting to note that the last two structures are gauge invariant 
provided we simultaneously contract them both with incoming and outgoing photons q\ v q2pX^ p = 
0. Namely, e.g. for the last tensor we get q% v T^) p = -\iV pv P v A^ and q 2fl T$ p = P p V pu Af, 
where we use throughout the convention A^ = ie pa A p and introduced a projector V pv = g pu — 
qx^q-ivlqx ■ q% which gives iP p V pv = q u + §(2f - r])P u and iV pv P v = q p + §(2£ + r])P p up to 
terms of order 0(A±) which we drop since they would exceed the accuracy we work to. A similar 
result we get for T^) with Aj_ being replaced by Aj_. Using the definitions of the photon and 
hadron momenta in terms of light-cone vectors, consequent simple algebra leads to the following 
structures in photon and hadron momenta 

^ / 

9pu 




Tffi — \te\poPp q<j 9fie 





qi ■ q2 



r it] P = i£e\ P *q* ( 9,,e - JT^J - ) , (18) 

where a missing twist-four part is restored minimally according to previous discussion. All struc- 
tures have a well defined forward limit. It will be shown in section [5] that these terms are equivalent 




to those used in our previous studies 0, p| . In subsequent sections we derive relations between 
the 'transverse' generalized distributions introduced here and conventional twist-two ones as well 
as three-particle correlation functions. 

3 Twist decomposition. 

In this section we present a decomposition of the two-quark operators into separate twist compo- 



nents to twist-three accuracy. Similar analyses have been done in Refs. W\ [13], 111] . However, an 
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essential new ingredient of our study is the treatment of operators with total derivatives. Since 
the group-theoretical notion of twist as the dimension minus spin of an operator is well defined for 
local operators, our strategy will be thus: first, an expansion of a light cone operator in infinite 
Taylor series in local ones, second, an extraction of definite twist components and as a final step 
the resummation of the result back into nonlocal form. An alternative approach directly based 



in terms of light-ray operators is described in Ref . [jn| . In the following we will not take care of 
trace terms, proportional to n p , since they only contribute at twist-four level. 
The Taylor expansion of Eq. ([11]) in terms of local operators simply reads 



O p (k, -k) 



3=0 J 



IK) 



-n pi . . . n^.Cp^j...^, 



with 



O 



p;p,i-Hj 



^TpiV^ ...iV^ip, (19) 



where V p = V p — V p and V p = d p — igB p . To have a well-defined decomposition into twist-two 
and three contributions we extract tensors 7Z p - pl ... P:j with (j + 1) indices corresponding to the two 

. This can easily be done with the result 



Young tables 



p 




1-1-2 




Pn 


and 


mi 


1-12 








p 





O p {k,-k) = K%k,-k)+K 3 p {k,-k) 

00 1 -ZK^ 



j=0 J' 



n P-i ■ ■ ■ n H j^Ui.-Ay 1 ■ + i'"p;ni... l x j I > 



(20) 



with 



K^...h = A S H $r p iv^ . . v>, (2i) 



where S is a symmetrization (and trace subtraction) operation of j indices with weight \ and 
antisymmetrization being defined by A t PlP2 = ~ (t piP2 — t P2Pl ). In terms of light-ray operators 

/^l — 

the twist-two part reads 



K 2 Jk, -«) 



du\ ip(—uK,n)r p i/j(Kn) 

dp (unn) — dp (—uku) — 2igB p {rnn) 



+ — j drip (— uku) R 



V>(u«n) \ (22) 



in the light-cone gauge B + = 0. If we would work in a covariant gauge this would result, apart 
from restoration of the gauge-link factors, to the substitution in the square brackets 



dp (uku)- dp (-unii) - 2igB i 



'p(rKn) — > V p {uku)— V p (—uku) + 2i tk gG +p (rKn). (23) 



One can easily project onto the twist-two part from the operator O p by contraction with the vector 
n p , namely, 1Z\ = + . 
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The conversion of the twist-three part in terms of three-particle operators is more tricky as 
compared to the forward case due to relevance of operators involving total derivatives. From the 
Eqs. (20) and (^) one can easily derive, however, the following equation, cf. 



V O p (k, -k) 



v TZ 2 p (k, —k) + k J q duliu e +ppv d p A O u [uK, —uk) 



(u — t) v S^(uk, tk, —uk) — (u + t) v S (uk, tk, —uk) 



(24) 



with d p = d P + d/i being the total derivative, and similar expression for A O p which is obtained 
upon substitution F p — * T p ^. These two equations set an infinite iteration in total derivatives. 
We can represent these relations in a schematic form 

v O = V 1Z 2 + J * (ed a O + V S) , a O = A TZ 2 + J * (td v O + A S) , 

with J being an integral operator. A formal solution to this system of equations reads 

v o = (i-j 2 d 2 y 1 *{ v n 2 + edj* A n 2 + j*(jd + i)* v s}, 

and similar for a O. So that l O — ! 1Z 2 is a twist-three part of i O. The series can be summed 
up into total translations in the light cone formalism. Let us demonstrate first the effect of total 
derivatives in the basis of local operators. The moments of Eq. ( P4"D look like 



Vq . = V K 2 



J 1 k=i 



(25) 



where obviously O 



p-,3 ~ n Pi 



v s± hk = i g v~ 2 i, dt l [i+g +p ± i 1+1 ,g +p \ ef-v 

and the analogous one for odd parity, i.e. V — > A and 7+ — > 7+75. The solution for 7Z 3 reads 
1 

Vn l;j = 0^ EO' - { a m ie +P^ ^ M AlZ l;j-l-i - °i ( id P n ° - 9 P Jd + ) v-Rl.^f A 

^ 1=0 

- 7 e'e W {0 - * - 0^+ - (-l)'*^^} , (27) 

J z=o fc=l 

with ct; = -[1 — (—1)']. Now we can resum this equation back into light-ray operators, or this can 
directly be obtained from Eq. (^5|), with the result 



n M .O p;Atl ... M . and we have introduced three-particle local operators 



Tfc-i 



(26) 



K 



+ K / G?T 



(«, (u - u)k) + A TZ 2 U ((u - u)k, -k) 
-w (<9 p ?v - # p(J <9 + ) y 7^ («, (u - u)k) - v K 2 a ((« - -«) 
(w — r) v Sp (k, (t + u)k, (u — u)k) — (u + t)^^ (( u ~~ u)k, (t — — k) 



(2* 



and same for V <-» A. Contrary to the forward scattering there appeared contributions from 
different parity operators to a given twist-three one and the 'center-of-mass' of two- and three- 
particle operators gets shifted by a total translation exp {±iuK&). This equation gives a relation 
between skewed parton distributions of different 'twists' when sandwiched between hadronic states, 
which we discuss in the next section. 



4 Twist-three skewed parton distributions. 

In this section we define twist-three skewed parton distributions and their relation to twist-two 
ones. Before we deal with spin-| functions, we consider first a more simple case of spinless target. 
For generality, we deal with the incoming and outgoing hadrons of different masses Pf = M\ ^ 
Pi = Ml 



4.1 Spin-0 target. 

It is instructive to start the analysis with the expectation values of local operators. Since for spin- 
hadron we can not form a 'twist-two' axial-vector, only the vector twist-two operator develops 
non-zero reduced matrix elements, which are given by 

(AlXm... J p i> = JS^ {/'„... P N B 3+W + A p l) K ■ ■ ■ P^Bj+u + • • • + A, . . . A H B j+lfi } , 
<AlX w .../J P i> = °" ( 29 ) 
The reduced matrix elements Bjk are defined as moments of a skewed parton distribution B(x, 77): 
1 d k t 1 



B 



3,3 -k 



k\ dt] k 



J dx x J 1 B(x, 77)1,7=0, where < k < j, 1 < j. (30) 



For the analysis of the moments of skewed parton distributions we have to project both sides of 
(|29|) with the light-cone vectors n Pl . . . . Since 



n ^ • • • n ^ J=L A " • • • A «*Vi • • • p »> = JTT { kA f> A + lp i +1 ~ l + + 1 - k)p p A k + pt k } 

= {kA^ 1 + + k)P pV k } Pi, (31) 

the coefficients in front of the two monomials are generated by derivatives of Bj + i(rj) w.r.t. the 
skewedness parameter, namely, 

m. d \ Pl ) = Pp Pi (1 - ^) b 3+m + \H^B j+ ^ 



{P2\ V K^\P l )n p = P{ +l B j+1 ( V ), (32) 
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with 

i+l 



B j+i(v) = J2 7 l k B j+i,j+i-k = dx x 3 B{x,r}), 

1 n J —1 



fc=0 

The moments Bj + i(r]) can be taken from the + component of the operator TZ + . As a consequence of 
symmetrization the first equation contains a Wandzura-Wilczek term proportional to A^ = A p — 
rjPp, which effectively enters as a twist-three contribution to the scattering amplitude. The matrix 
element of the light-ray operator can be obtained in a straightforward manner by resummation: 

(P 2 \ v H 2 p {k,-k)\P 1 ) = j^dxe- iKP + x (p p B(x, V ) + Aj J^dyW 2 (x,y)^-B(x, V ) ) j , (33) 

where the kernel reads W 2 (x, y) = 8(x)6(y — x)/y + 9(— x)9(x —y)/ (—%)■ 

Next we define the reduced matrix elements of the antiquark-gluon-quark operators. Since 
these operators are partially antisymmetrized, we have obviously two vectors A^ and the dual 
ones Ap- = ie^A^, introduced above, in our disposal. Thus, the general decomposition of reduced 
matrix elements reads 

(P*\ V S%, k \Pi) = K S tkPl WS% >k \Pi) = A;/;,,/". (34) 



The duality relation (|T3|) reduces immediately the number of independent contributions to two 
instead of four, namely, Rf k = Sj k . It turns out convenient to work in a mixed representation 
for the skewed parton distributions. We introduce a representation that depends on the position 
of the gluon field and a Fourier conjugate variable with respect to a distance between both quark 
fields: 

(p*\llt[ K : UK :~i\\pi) = pi\ f p \ i <i.«- s-u,,,,). (30) 



4S^(/c, uk, — k) j |A 



The moments with respect to the momentum fraction x are given by a polynomial of order j — 2 
in the variable u: 

?(«^) = //«^(x,« j)? ) = E(f|) {^y^^T^ S tk(v), (36) 

where Sj t k(i]) are polynomials in rj of order j defined in Eq. (34). As a simple consequence of our 
definition we have 

fju^Sfiu^) = 2£j^£fa{v), jju^SfM = (37) 
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Finally, to find the matrix element of the operator O p j 
findings ( |3~2"1) and (|37|) into the solution (|27|) : 



remains to insert our 



A 



P?{P 2 \ v O p , d \P x ) = [P p + -f 2J i+1 fo) J 



fc=0 



7 + 1 



drj 



. v y-k [ — -—7— ] Bm-iW 



77 



fc!(-7?) 



" P ^ 2 (j + l)(^-2)!7-i 



l —^st{u )V ) - {-iy- kl -^s- k {u^) \ . 



(3* 



Now we consider the axial-vector case, which is handled in the same way. The main difference 
is that there is no twist-two part for spinless target. However, a non- vanishing twist-two part of 
the vector operator induces now a Wandzura-Wilczek type relation due to the epsilon tensor in 
Eq. ( p7|) and provides 



p: 



' cfj-k , ,4-h ( d k + 1 



drj i] 



(39) 



k\(-rj) 



k=2 



1 — It 



— Ajy - ^ y ''' / du\^^Rt(u, V ) - (-l) j - k] ^RZ(u,r)) 

p ^ (j + l)(fc-2)! 7-1 12 fcV /; v ; 2 " ' 



The final step is a summation of local operators, see Eq. ( pOf) , which leads to the expectation 
values of the light-ray operators in terms of Fourier transform of skewed parton distributions in 
parity even and odd cases 



-■>r,,, Ai l l p. w JZV\t*_V* ]B(Si , ) 



— 1 [77 1 V 77 T) J \ drj 77 dy 



(40) 



+ \ [ P p + -A p L jB(x, V ) 



2 V 7 / vj 



2 \ V VJ 



{P 2 \ a O p {k,-k)\P x ) = I" dxe^Uj f %V-[-A (T---T-I *<•"•'/> 

J-i I p J-i \r]\ \rj rj J \ drj rj dy 



(41) 



1 |r/| J-i 



2 U 77/ 2 V 77 77/ 



Here W" f ± £ , ± a l = ■f s W(±^, ±*) and the W kernels read 
V. v v J d y \ v' vJ 



w(x,y) = ^^, with o( Xi y) = S i gn (i + y)e(^Ae^ y x 



l + y 

W±(x, y) = l - {W(x, y) ± W(-x, -y)} . 



l + y \l + y r 



(42) 
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In restoration of non-local form one uses the following result for Mellin moments of ^-kernels 

e„ i + i M U V fr„ 20 + 1) 







- 













fc=o J 1 ■*• I'/l \v '// fc =0 

(43) 

Note that the piece ^A p l B(x, rj) in Eq. (|40| ) is cancelled by a term arising from the convolution 
with the W + kernel. 

4.2 Spin-i target. 

Now we are in a position to discuss a spin-- target. It is convenient to express the expectation 
value of local operators in terms of spinor bilinears 

(b,b) = U{P 2 ,S 2 )(l, l5 )U(P 1 ,S 1 ), (h p ~h p ) = U(P 2 ,S 2 ) lp (1, 75 ) U(P l ,S 1 ), 
(t P aJ P a) = U(P 2 ,S 2 )ia p<7 (l, l5 )U(P 1 ,S 1 ), (44) 

Obviously, the dual tensor form factor t pa is obtained from t pa by contraction with the e-tensor 
and can, therefore, be eliminated. Furthermore, equation of motion shows that in each parity 
sector we have the relation between the structures (0) 

P p b = M + h p -t pa A°, A p b = M_h p -t pa P°, (45) 
A p b = M + h p - t pa P\ P p b = M_h p - t pa A\ 

where M± = M 2 ± Mi. As demonstrated above for a scalar target, the symmetrization provides 
us a Wandzura-Wilczek term proportional to A -1 . To take advantage of the analysis already 
performed in the preceding section and for a symmetrical handling of even and odd parity sectors, 
our basis is spanned by the (pseudo) scalar and the (axial) vector bilinears. At the end, we express 
the result in terms of conventional ones, introduced by Ji ||, 

A CT ~ A ~ 

hp, e p = tpa — , h p , e p = j-^-b. (46) 

Obviously, for the (pseudo) scalar bilinears we can take the results deduced for scalar target. The 
only new structure for the spin-| target is proportional to the h p , h p . 



In parallel to section |4.1| , the local matrix elements read 



(P2\ V K, H ..J P i) = „S W h„ {P m . . . P„ t A j+1J+1 + . . . + A w . . . A n A j+lfi } 

where the ellipses (here and later on) stand for the r.h.s. of Eq. (p9|) (and corresponding equations 
from the preceding subsection). Analogous relation holds for the parity odd case. The moments 
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of the skewed parton distribution are defined analogously to Eq. (|30|), Aj + i(rj) = j\ dx x-M(x, r]). 
The difference arises only from the symmetrization procedure 

n w . . . n w S ^A M1 . . . A Mfc P Mfc+1 . . . P = Wp + + (j - k)P p h + rj k + fcApM* -1 } ,(48) 
and provides now a new structure 

(^lX;l p i> = j^j (tf + + - + K K-^j A j+1 (v) + ^ • • • , 
(AlX^iK = ^i+ifo) + B i+M- (49) 



Next we resum the local result and get 



(P 2 \ V K% K , -k) {P.) = £ dxe- iKP Ap p (j±A(x, rj) + ^Bfo rj) 



(50) 



i 

dy W 2 (a;,2/) 

i 



Projecting this expression with vector n p we obtain conventional definitions for twist-two skewed 
parton distributions. The basis used presently can be easily expressed in terms of Ji's parametriza- 
tion as follows, for Dirac structures 

6=-^+, b = ^(h + -e + ), (51) 

and skewed parton distributions 

A = H, B = r)E, A = H + E, B = -E. (52) 

Po 1„ \ (_L Po 



For the following, it is useful to note as well that M + [h p — p^/i+J = \ t pa — ^-t +(T j A CT . 

For the time being we introduce three-particle skewed parton distributions without an explicit 
spinor bilinear decomposition, 

A discussion of the parametrization will be given below. 

Due to length of the consequent formulas we give below only results for the vector case since the 
axial one follows from the former by substitutions. Combining the Eq. P5|) with the two-particle 
part of Eq. fl2"T| ) we get for the /i-part of the local operators v O p 

Pl->WW = (P, + f) h + A j+1 + ± ^(-,)^{a >+ [± - k -±l) ,54, 



k=0 



+ iej a (h a P+ - P a h + ) jA fc+ i + 
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Transforming it finally to the non-local form, and adding the 6-term and the missing three-particle 
piece, one finds 



[P 2 \ v O p {-k,k)\P x ) = j\xe- iKP + x { (p p + ±A^ D( 



x, rj) (55) 




W" I * I fij(y, it, 77) + 1±V" (--, -1L] S;(y, it, ^ 
where we have introduced a shorthand notation for the combinations 

D (x, V) = jr A (%i v) + J^ B ( X ' V), D (Vi V) = jr A (y, rj) + jj-B(y, rj). (56) 

The derivatives w.r.t. rj in Eq. (|55"D acts only on skewed parton distributions but not on spinor 
bilinears. The axial-vector case is deduced by the following trivial substitutions: b b,h p <-> 
/t p , 5 <-> 5, A «-> A and 5^ — > i?^. Conventional form is obtained by means of substitutions 
( |5T| . [52| ) . Equation (^) is our final result. Its transverse part gives an expression for the 'transverse' 
twist-three skewed functions in terms of 'known' leading twist and three particle functions. 

Let us address briefly the parametrization of genuine twist-three part of the correlation func- 
tions. An analysis shows that there exist four independent twist-three spinor bilinears, which we 
can choose to be 

A^ — A^- A 1 — A^- (57) 
A "M + ' A 'P + ' A 'M + ' A 'P + ' 

while the dual ones read in the same sequence 



A^-L ft-L^f. A x — A^ (& 
A 'M' A 'P' p Mi' ^P/ [bt 



Note that a fifth candidate, M+j^, and trace terms proportional to M^_n p enter at twist-four 
level. One can immediately see that bilinears used in Eq. (|5"5j) can be spanned by this basis via 
relations h p - P p ^- = - AA^ + i £f |± and L - P p ^± = --^A^ + * A^. The 

h P F-f. 1— r?- 3 P Ff 1 — 77^ P F+ v f F+ 1— r/- 3 P F+ 1— rj-* P Ff 

parametrization of correlation functions S and R then reads 



Rf(x } u, V ) = A^Rf + Aj^Rt + A^Rf + Aj^Rt (59) 
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where duality (|13|) requires = Rf, Rf = Sf. To reexpress the parity even and odd sector only 
in terms of vector and axial-vector form vectors, respectively, one can equally use an alternative 
basis, where the independent elements are 

Af— , h p -pj^, Af— , L-pkh, (60) 
p M+ p P P+ p M+ p P P+ V ; 

while the dual ones read in the same sequence 



5 Power suppression of twist-three effects in DVCS. 

Now we point out the phenomenological consequences of our analysis for the DVCS process, 
eN — > e'N'j, with 77 = — £. Since DVCS interferes with the Bethe-Heitler process, the possibility 
exists to get a direct access to the skewed parton distributions via the interference term. In leading 
twist-two approximation of the DVCS hadronic tensor the amplitude squares for unpolarized or 
longitudinal polarized nucleon behave as: 

l^vcsl 2 oc ij, T BH T* VCS cx -jJ^ (l - ^f) ' , |T BH | 2 cc -L, (62) 

where A^ in = — 4M 2 £ 2 /(1 — £ 2 ) follows from the kinematical boundary on which A 1 - vanish. 
Consequently, the interference terms and so also the charge and spin asymmetries vanish (this is 
not the case for transversely polarized target). Since the leading Aj_ dependence could also arise 
from the twist-three contributions to the DVCS amplitude, one may argue that those enter in the 
interference term without power suppression. Consequently, the whole twist-two analyses would 
be spoiled. This possible complication has been already studied in the past [IT], unfortunately, it 



was not clearly emphasized that twist-three contributions are power suppressed. 

Our aim is to show that this is indeed the case, while the complete interference term will 
be published elsewhere |T3|. In our analysis |TTJ we used the following parametrization of the 
hadronic scattering amplitude that arises from the operator product expansion in the free fermion 
theory, where current conservation was restored in calculations by means of projection operators: 

T^(q, P, A) = -v g^v^t-L^L + {V^P a V pv + VppPvVn,) (63) 

p . q P ■ q z 

-<p i e <p Alf. --p i e p p ^p. 

I UXT (TTQp ' 7T n ' Ucr *- (TT F p ' TV j-, ) 

P ■ q P ■ q 

where has been introduced in section |2|. Comparing with the results given in Eqs. (|!6|]l8l) , we 
find that the form factors are related to each other in the following manner 

V llx = J jC^(x,Z) V 0»(x,ri), A lfi = J jC^(x,0 A O,(x,r]), A 2 , = 0, (64) 
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Vo 



1 Pn 



q- v O(x, V ))+C^\x,0 



1 £/j,vAq A 

2P-q 



O u (x,v) 



2P -q 

A straightforward calculation and power counting show that in the case of DVCS kinematics, 
where r\ = — £, the contributions of Vi and A 2 are power suppressed in comparison to the known 
twist-two result: 



^BH^DVCS 



2~2y + y 2 g 
1 - y A V 

1-y AVI 



k n 



y , 



J.- a, 



q ■ V{ + ie kqAJ 



•4 



k a - -q c 

y 



(65) 



Here J M is the electromagnetic current, A and denote the polarization and momentum of the 
incoming electron. We employed a simple rule that g M , k^ — ^q^, and P^, A^, when contracted with 
an electromagnetic current or spinor bilinears, give terms of order Q 2 , Q 1 , and Q°, respectively. 
Note that in agreement with this counting our results for unpolarized nucleon target coincide with 



Ref. ||16|| , where the contribution of V2 has been taken into account. Let us remark, that also in 
the case of a scalar target the same situation holds true, as it already has been shown in 0. 

Since, the twist-three contributions enter as 1/Q power suppressed term to all single spin 
and charge asymmetries, the comparison of models for twist- 2 skewed parton distributions with 
experimental data will be contaminated in an expected manner. Moreover, we emphasize again 
that the twist-3 functions are completely known in terms of twist-two ones in absence of the gluonic 
contributions as it is the case in the naive parton model. It is an interesting non-perturbative 
problem to estimate the size of the gluonic contributions in comparison to the Wandzura-Wilczek 
term. Assuming that the dynamical twist-three effects encoded in the three-particle operators are 
small, which presumably happens in view of recent experimental data |T7j and new lattice results 
18 1 for the forward kinematics, one can use the Wandzura-Wilczek part of the relation (|55|) as a 



model for the 'transverse' twist-three skewed parton distributions, {P2\ i O^\Pi) . This will allow to 
give a numerical estimate of the power suppressed contributions to asymmetries for the kinematics 
of present experiments Jl5 . 



6 Logarithmic scaling violation. 

Let us add a final remark on the logarithmic scaling violation of twist-three functions. Obvi- 
ously the Wandzura-Wilczek part evolves via the familiar twist-two generalized exclusive evolution 
equation known to two-loop accuracy nowadays [19]. The three-particle piece does not require a 



new study as well, since it was elaborated in great detail recently in the context of twist-three 
functions measured in inclusive reactions (forward limit restrictions was not made in those stud- 
ies). Namely, the operators S alluded to above fall into a class of the so-called quasi-partonic 
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ones |2(J and thus at leading order in the coupling constant their evolution kernel is given by a 



sum of conventional skewed twist-two ones with appropriate quantum numbers in subchannels: 



Kqg q = Kqg + K gq + Kq q . The issue of its diagonalization (see a recent review |21]) has been 
addressed in detail recently |23], |23J] within the context of integrable open spin chain models 
which arise in multicolour limitQ Namely, for N c — > oo the kernel K qq can be neglected and, 

e.g. jj-K-qg = V(4 + D + ^(4 - f ) - 2^(1) and = </>(4 + §) + V>(4 " |) ~ 2 ^(1) 

for <S + , and with | and | being interchanged for S~ . Due to conformal invariance at tree level, 
the kernels depend on the quadratic Casimir operator of the collinear conformal group SO (2, 1) 
or, in other words, on the conformal spin in the subchannel only, J = J(J — 1). It turns out 
that antiquark-gluon-quark system admits an extra integral of motion Q$ and is thus completely 
integrable. The result of the analysis allows to find the eigenfunctions & and eigenvalues S of the 
system and thus solve the evolution equation, 

(a s (Ql)\ NcS{a}/P ° 3 
S(x 1 ,x 2 ,x 3 \Q 2 ) = ^2^ {a y(x 1 ,x 2 ,x 3 ) s ((S {a} (Ql))), 22xi=r). 

{ a } \as{Q 2 )J i=1 

Here /3o = ^TpN f — g-C^ is the QCD beta function and a set of quantum numbers {a} parametrizes 
solutions and can be chosen as eigenvalues of the conformal spin J and the charge Qs- Namely, 
the lowest anomalous dimension, Sq(J), is known exactly [p5[] , while the rest of the spectrum, 



S q (J), can be described with a high accuracy using WKB approximation [22, 23. 2_4] with the 
results 

£ (J) = ij)(J + 3) + *0(J + 4) - 2^(1) - \, £ Q {J) = 2 In J - #(1) + 2 Re if, + i Vs \ - | 



where r/s is related to the conserved charge by the relation rjs = \\J2Q S / J 2 — 3 and obeys a WKB 
quantization condition which, once solved, gives quantized values of Sq. 

For practical purposes, one may also generalize the polynomial reconstruction method as it 
was applied for the forward kinematics in Ref. [^6] or use a brute force numerical integration to 
solve the evolution equations. 



7 Conclusions. 

In this paper we have studied twist-three effects for the two-photon amplitudes in the generalized 
Bjorken kinematics. We have demonstrated explicitly the restoration of the electromagnetic gauge 
invariance for a Compton-type amplitude to twist-four accuracy. We have given as well the exact 

3 The latter gives a good approximation since the corrections to the leading result are suppressed in l/N^. 
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Lorentz structure of the amplitude motivated by our results at leading order in the coupling con- 
stant. Obviously, beyond leading order, e.g. different Lorentz structures in T^) will be multiplied 



by independent functions (see Ref. [10]), analogues of F\ and F2 in the forward scattering. Our 
analysis demonstrates a necessity of introduction of twist-three two-particle generalized distribu- 
tions. The latter can be related by means of QCD equation of motion and Lorentz invariance to 
the familiar twist-two ones and interaction dependent antiquark-gluon-quark correlation functions. 
The former ones represent a generalization of the Wandzura-Wilczek type relation, while the last 
ones give dynamical twist-three contributions. With the assumption of smallness of dynamical 
contributions the Wandzura-Wilczek part with parametrization of leading twist skewed parton 
distributions can serve as a model for 'transverse' functions. Contributions of these functions to 



diverse asymmetries [1 1J will be considered elsewhere |T5|. 



Note added: Recently there appeared a note |27| where the quark hand-bag diagram is 
calculated in parton model with transverse momentum being kept. This result overlaps with a 
part of our analysis in section [2|. 

We thank A. Kirchner for a collaboration at an early stages of the work. A.B. would like 
to thank A. Schafer for the hospitality extended to him at the Institut fur Theoretische Physik, 
Universitat Regensburg. This work was supported by Alexander von Humboldt Foundation, in 
part by National Science Foundation, under grant PHY9722101 (A.B.) and by DFG and BMBF 
(D.M.). 
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